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ABSTRACT 


-^Approximately  2000  kilometers  ( ~  1250  statute  miles) 
of  re f lectome ter  data  collected  within  160  kilometers 
(100  statute  miles)  of  the  ice  edge  in  the  North 
American  Arctic  were  analyzed.  The  ref lectometer  sig¬ 
nal,  which  shows  a  sharp  decrease  in  areas  of  open 
water/thin  ice,  was  used  to  initiate  and  develop  a 
method  to  begin  an  evaluation  of  the  frequency  of 
occurrence  and  percentage  of  open  water  from  the  ice 
edge  to  approximately  160  kilometers  (100  statute 
miles).  Comparisons  were  made  within  and  among 
regional  data  sets.  Individual  regions  were  not  unam¬ 
biguously  identifiable  by  lead  width  and  frequency 
characteristics.  Distance  into  the  pack  from  the  ice 
edge  did  not  have  a  direct  relationship  to  the  fre¬ 
quency  or  percentage  of  open  water.  The  result  of  no 
apparent  relationship  between  the  frequency  of  occur¬ 
rence  and  percent  of  area  of  open  water  may  be  due  to 
the  restricted  samples — restricted  in  season  and 
total  area  covered. 
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OPEN  WATER  AND  THIN  ICE  DETECTION  IN  THE 


ARCTIC  MARGINAL  ICE  ZONE 
USING  REFLECTOMETER  SIGNAL  ANALYSIS 


INTRODUCTION 

The  extent  of  Arctic  ice  coveraje  ranges  from  a  summer 
average  minimum  of  5.2  million  kin^  to  a  winter  average  maxi¬ 
mum  of  11.7  million  km2  (wittman,  1959).  The  precise  extent 
of  the  ice  edge  at  any  time  is  a  result  of  previous  meteoro¬ 
logical  and  oceanographic  conditions.  In  the  course  of  its 
southward  growtn  ana  drift  during  winter  montns  tne  outer  mar¬ 
gin  of  the  Arctic  pack  eventually  merges  with  sheets  of  new 
fast  ice  growing  seawaro  from  tne  coastlines,  ending  surface 
navigation  for  all  ships  except  icebreaker  assisted  operation. 
These  operations  are  primarily  along  tne  fringe  of.  the  pacx 
with  limitation  of  mobility  increasing  with  pack  penetration. 
From  October  to  June  the  Arctic  Ocean  remains  virtually  ice- 
locked;  however,  tides,  winds,  and  currents  can  produce  areas 
of  open  water  within  tne  pack  and  along  the  marginal  ice  zone 
at  any  time. 

The  purpose  of  this  investigation  is  to  develop  a  metnoo 
and  to  oegin  an  evaluation  of  the  frequency  and  percentage  of 
open  water  from  the  ice  edge  to  approximately  100  miles  into 
the  pack. 

Open  water  associated  with  the  ice  edge  can  be  placed  in 
one  of  five  categories.  A  crack  is  a  small  unnavigable  break 
caused  by  tides,  temperature  change,  current,  or  wind.  A  leau 
is  a  long  narrow  navigable  water  passage  in  pack  ice  and  a 
polynya  Is  any  sizable  sea  water  area,  other  than  a  leau, 
encompassed  by  sea  ice.  A  bay  is  a  minor,  wi.tn  a  bight  a  major 
inward  bend  of  the  ice  edge  or  ice  limit  formed  either  by  wind 
or  cjrrent.  Because  of  the  fine  resolution  obtainable  in  tiiis 
study  (15  meters),  a  crack,  lead,  and  polynya  will  De  con¬ 
tained  in  one  category  with  a  Day  or  bight  being  regarded  as 
the  ice  edge.  Average  ice  limits  were  plotted  relating  to  tne 
year  and  season  of  data  collection  and  areas  of  open  water 
were  analyzed  from  tne  ice  edge  to  approximately  100  miles 
into  the  pack.  Seasonal  ice  limits  have  been  estimated  Dut  for 
increased  precision  tne  reader  is  directed  to  tne  Eastern- 
Western  Arctic  Sea  Ice  Analysis  annually  prepared  oy  tne  wavai 
Polar  Oceanograpny  Center  (NPQC),  Suitlano,  Md .  ,  wnicn  was 
used  in  this  study. 

Satellite  data  interpretation  over  the  Arctic  regions 
shows  major  areas  of  open  water,  but  the  low  resolution  limits 
spatial  detection.  Thus,  precise  percentages  of  water  witnin 
the  pack  are  difficult  to  obtain  by  this  method. 
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The  two  types  oC  data  combined  for  this  study  are  aerial 
photography  and  ref lectometer  signal  analysis  in  conjunction 
with  a  laser  prof  i  iometer .  The  data  has  been  collected  over 
the  Arctic  pack  ice  on  an  opportunity  basis  since  1970  using 
RC-8  and  RC-10  aerial  cameras  and  a  Spectra  Physics  Geodolite 
3A  laser  terrain  profiling  system,  which  uses  a  modulated  con¬ 
tinuous  wave  laser  technique  to  obtain  a  precise  measurement 
of  instrument  height  above  the  surface.  Details  of  laser 
analysis  and  a  description  of  the  system  can  be  found  in 
Ketchum  (1971)  and  Welsh  and  Tucker  (1971). 

The  laser  signal  is  stored  on  magnetic  tapes  and  is  used 
to  classify  dynamic  ice  parameters  such  as  surface  roughness, 
ridge  height  distribution  and  frequency,  and  power  spectral 
density,  as  well  as  for  input  to  and  refinement  of  various  ice 
prediction  models  of  the  Arctic. 

Three  other  signals  are  recorded  on  coincident  channels 
to  compliment  laser  data,  which  are  time  code,  phase  lock 
fail,  and  ref  lectometer .  When  the  time  code  record  is  corre¬ 
lated  with  the  aircraft  navigation  logs,  accurate  track  lines 
are  reconstructed.  The  re f lectometer  channel  records  the  meas¬ 
ured  millivolt  change  in  light  intensity,  which  is  the  total 
of  the  laser  light  and  the  sunlight  reflection  passing  through 
a  3  angstrom  optical  filter  centered  on  6328  angstroms.  Open 
water  is  seen  as  a  sharp  decrease  in  the  reflected  laser  light 
intensity  (Fig.  1)  (Wilheit,  Nordberg,  and  others,  1972).  This 
figure  will  be  discussed  in  detail  in  the  later  section  on 
data  analysis.  The  ref lectometer  signal  is  "noncalibrated; " 
therefore,  all  measurements  are  relative  to  surrounding 
values.  The  phase  lock  fail  channel  indicates  loss  of  laser 
signal  caused  by  envi ronmental  conditions,  i.e.,  clouds.  This 
signal  can  be  used  as  a  check  to  insure  that  the  laser  record 
has  not  been  geographically  distorted  by  environmental  condi¬ 
tions.  A  detailed  description  of  the  laser  data  reduction 
process  written  and  utilized  at  NORDA  (Naval  Ocean  Research 
and  Development  Activity)  can  be  found  in  Lohanick  (1981). 

DATA  ANALYSIS 


The  selected  magnetic  analog  tape  containing  the  chosen 
data  track  is  played  on  an  Ampex  FY-1300  tape  recorder  through 
an  HP2240A  (analog/digi tal  converter).  The  analog  voltages  are 
then  digitized  to  make  them  compatible  for  reading  by  the 
HP984513  (tabletop  computer).  The  data  is  stored  on  flexible 
discs  driven  by  an  UP9885M  (flexible  disc  drive)  for  manipula¬ 
tion,  and  a  plot  of  the  data  is  generated  for  analysis. 

The  preliminary  stage  of  this  investigation  involved 
matching  aerial  photography  with  the  plots  of  the  ref lectom¬ 
eter  signal  to  determine  the  relative  change  of  signal  inten¬ 
sity  with  ice  thickness.  \  distinctive  drop  in  signal  strenjth 
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occurred  in  areas  of  open  water  and  thin  ice.  Thin  ice,  for 
the  purpose  of  this  study  is  defined  as  <30  cm.  As  shown  in 
Figure  1,  the  change  and  variation  in  signal  intensity  is  a 
direct  measurement  of  surface  albedo.  In  newly  forming  ice, 
the  surface  albedo  will  generally  have  a  direct  correlation 
with  thickness.  In  Figure  1,  the  ref lectometer  transition  at 
point  A  from  a  homogeneous  lower  return  to  a  higher  response 
with  increased  variation  is  shown  on  the  photograph  as  a  newly 
ridged  zone  from  A  to  B .  The  signal  change  at  point  C  is  not 
as  great  or  distinctive  as  at  point  D  because  the  newly  form¬ 
ing  ice  evident  in  the  photograph  increases  the  albedo  and  is 
starting  to  blend  the  signal  to  surrounding  features.  The 
distinction  at  the  lead  boundaries  (points  C  and  D)  allow  for 
the  fine  resolution  of  this  procedure.  Point  E  shows  the 
signal  response  across  a  ridge;  the  great  variation  in  surface 
reflectivity  due  to  roughness  may  lead  to  another  future 
method  of  discerning  ice  types  with  an  optical  system.  These 
data  plots  are  scaled  according  to  the  aircraft  navigation 
logs,  and  the  width  and  frequency  of  these  areas  are  catalog¬ 
ed  for  analysis  (Appendix  B).  Computer  programs  were  written 
at  NORDA  for  this  investigation  to  automatically  identify 
areas  of  open  water/thin  ice  in  the  re f lectometer  signal  but 
the  methods  proved  inadequate  or  unworkable.  Holyer  et  al. 
(1977)  discusses  some  problems  of  automatic  data  analysis  with 
the  laser  signal.  These,  along  with  a  "relative"  rather  than 
"calibrated"  signal,  complicates  the  procedure;  therefore, 
more  effort  will  be  required  in  the  future  to  generate  a 
reliable,  totally  automated  procedure. 

STATISTICAL  PROCEDURE 

The  raw  data  from  the  ref  lectometer  signal  yields  the 
frequency  and  width  of  areas  of  open  water/thin  ice  along  the 
aircraft  track.  A  Wilcoxson's  Sum  of  Ranks  Test  was  then  used 
to  determine  the  level  of  significance  between  and  within  the 
data  sets.  This  statistical  test  was  used  because  it  allows 
nonparametric  comparison  of  two  populations  based  on  indepen¬ 
dent  random  samples  and  is  insensitive  to  the  dispersion  of 
measurements  in  the  sample.  It  is  also  free  of  the  invalid 
assumptions  of  normality.  Testing  for  a  null  hypothesis  (H0) 
of  no  difference  between  the  samples  with  the  alternative 
(H^)  realizing  a  difference  to  a  measurable  level,  a  proba¬ 
bility  (P)  >  5%  is  interpreted  as  no  significant  difference 
being  proven  by  this  test,  P  =  5%  or  less  is  regarded  probably 
significant,  and  P  =  1%  or  less  statistically  significant. 
Original  data  rather  than  the  "class  interval"  data  shown  in 
Appendix  B  was  used  to  eliminate  tied  ranks  that  tend  to 
weaken  the  power  of  this  test.  A  detailed  description  of  the 
test  can  be  found  in  Langley  (1970).  The  statistical  test  was 
run  on  the  data  sets  using  a  HP9845B  (tabletop  computer).  A 
listing  of  the  program  written  at  NORDA  for  this  investigation 
can  be  found  in  Appendix  A.  Appendix  D  shows  a  sample  printout 
from  this  program. 
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Figure  I .  Relative  change  in  response  of  reflectometer  signal  as  a  lead  is  crossed 
showing  the  direct  relation  of  signal  intensity  to  surface  albedo 
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The  data  sets  were  separated  into  consecutive  5  minute 
segments  to  determine  levels  of  significance  within  regions 
and  were  combined  geographically  to  test  significance  between 
regions . 

RESULTS 

The  ice  conditions  described  in  the  following  regional 
summaries  are  sometimes  referred  to  as  "minimum"  or  "maximum", 
which  designates  the  extent  of  southern  growth  and  accumula¬ 
tion  with  no  rejard  to  ice  type  distribution.  Therefore,  at 
maximum  conditions  the  ice  is  generally  thinner  in  the  first 
100  miles  from  the  ice  edge  due  to  its  recent  formation  and 
thicker  during  minimum  extent  due  to  the  higher  percentage  of 
multi-year  ice  contained  in  the  boundary  zone .  The  following 
summaries  are  related  to  the  divisions  shown  in  Figure  2.  The 
laser/ref lectometer  data  was  collected  on  an  opportunity  ba¬ 
sis;  therefore,  seasonal  comparisons  are  limited.  The  dates  of 
data  collection  are  included  with  the  histograms  in  appendix  C 
and  will  be  seasonally  correlated  with  future  data. 

LINCOLN  SF.A 

Approximately  380  km  of  re f lectometer  data  was  collected 
in  the  Lincoln  Sea  on  6  November  1970  (Fig.  2,  track  line  5). 
The  data  was  analyzed  in  14  consecutive  5-minute  segments 
(approximately  27  km  per  segment)  originating  160  km  from  the 
ice  edge.  Ninety-one  2-digit  combinations  exist  with  14  data 
sets  but  only  86  combinations  were  tested  because  five  combi¬ 
nations  did  not  combine  to  form  at  least  10  elements,  which  is 
required  for  the  test. 

The  first  segment  analyzed  (segment  farthest  from  ice 
edge)  was  one  of  two  segments  found  to  be  significantly  dif¬ 
ferent  f c  xn  others  in  this  region.  The  probabilities  that  the 
data  from  this  segment  came  from  the  same  population  as  the 
other  segments  in  this  regions  were:  <0.2%  for  two  segments, 
0.2-1%  for  two  segments,  1-5%  for  4  segments,  5-10%  for  two 
segments,  and  >10%  for  two  segments.  This  shows  that  for  this 
data  set  a  change  in  lead  frequency  and  width  characteristics 
occur  about  100  miles  from  the  ice  edge.  Four  of  the  five 
segments  of  no  significant  difference  (P  >  5%)  for  this  seg¬ 
ment  occur  50-90  km  from  the  ice  edge,  which  identifies  a 
wi thin-region  variation  separate  from  the  total  population. 
Segment  10,  approximately  40  km  from  the  ice  edge  showed  a 
probable  significant  difference  from  two  segments  (P  =  1-5%) 
and  statistically  significant  difference  (P  <  1%)  for  two 
other  segments  with  no  apparent  relationship  to  pack 
penetration. 

In  summary,  74  of  the  86  pairs  of  data  sets  compared 
showed  no  significant  difference  with  respect  to  lead  width 
and  frequency.  The  first  100  miles  from  the  ice  edge  is  a 
relatively  homogeneous  zone  with  a  significant  boundary 
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occurring  beyond  this  distance.  There  is  also  no  apparent  re¬ 
lationship  between  the  percentage  of  open  water  witn  distance 
from  the  ice  edge  as  shown  in  Appendix  B.  Tne  region  nad  an 
overall  open  water  coverage  of  3.78%. 

NORTH  SLOPE  REGION 

The  North  Slope  region  contains  tne  Beaufort  and  CnuKcni 
Seas  and  is  one  of  the  most  studied  as  well  as  most  strategic 
areas  in  the  Arctic.  Three  track  lines  of  ref lectometer  data 
were  analyzed  as  shown  in  Figure  2.  Track  2  was  flown  13  No¬ 
vember  1970,  perpendicular  to  the  ice  edge  150  km.  Each  of  tne 
five  data  segments  were  compared  with  each  of  the  other  seg¬ 
ments  (total  of  10  comparisons),  and  only  one  pair  snowed 
statistical  difference.  This  pair  was  the  first  and  last 
segment  of  the  track.  Possibly  a  subtle  change  is  occurring 
with  each  segment  as  the  pack  is  penetrated,  and  it  requires 
100  miles  of  separation  before  the  difference  can  be  statisti¬ 
cally  observed.  This  track  line  crossed  the  ice  edge  (first 
segment),  which  is  an  unstable  area  and  had  10.99%  open  water 
compared  to  the  lowest  segment  of  1.11%  water.  The  first 
segment  beyond  the  160  km  limit  of  this  investigation  was  also 
statistically  different  from  the  segment  containing  the  ice 
edge  and  had  5.32%  open  water.  Therefore,  a  general  statement 
that  percentage  of  open  water  decreases  with  distance  from  ice 
edge  is  not  valid  in  this  population. 

The  second  track  line  was  flown  13  May  1971,  parallel  to 
the  coast  of  Banks  Island  130  km  from  shore  for  a  distance  of 
86  km  (Fig.  2,  Track  4).  This  area  is  one  of  the  rougnest 
areas  in  the  Arctic  due  to  ice  movement  in  tne  Beaufort  Gyre. 
The  area  had  mostly  small  fractures  (15-90  m),  and  the  data 
segments  averaged  only  2.74%  open  water  witn  no  significant 
difference  found  between  any  segment  comparisons. 

The  third  track  line  collected  data  off  the  Nortn  slope 
of  Alaska  during  maximum  ice  conditions  (24  April  1974).  Tnis 
track  line  paralleled  the  Coast  from  145°  to  156°N  longitude 
at  both  65  and  93  km  (Fig.  2,  Track  3).  Tne  individual  oat  a 
segments  are  not  contained  in  Appendix  B  because  no  water 
openings  15  m  were  found.  A  few  small  cracks  (5-10  m)  ap¬ 
peared  occasionally  in  the  data  with  no  apparent  pattern.  Tn<_ 
lack  of  open  water  in  this  data  set  is  a  result  of  a  general 
southward  drift  of  sea  ice  under  the  influence  of  tne  pre¬ 
vailing  northerly  winds  present  in  the  Beaufort  Sea  during 
this  time  of  year.  The  only  persistent  open  water  in  tins 
region  during  maximum  ice  conditions  determined  from  satellite 
imagery  is  off  the  south-facing  coast  east  of  Point  Hope. 

Wadhams  and  Horne  (1980)  analyzed  submarine  sonar  data 
collected  in  April  1976  in  the  same  area  as  track  3.  In  tneir 
study,  a  lead  was  defined  as  a  continuous  sequence  of  dept.) 
points  in  which  no  point  exceeds  1  m  in  draft.  Their  results 
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showed  that  98%  of  the  leads  were  <50  m  cross-track  and,  if  a 
submarine  requires  a  200  m  lead  for  a  safe  surfacing,  it  would 
have  to  travel  68  km  to  find  one. 

BERING  SEA 

Ref lectometer  data  in  this  region  was  collected  on  9  Feb¬ 
ruary  1976  (Fig.  2,  Track  1).  Beginning  160  km  from  the  ice 
edge,  six  5-minute  consecutive  segments  were  obtained  perpen¬ 
dicular  to  the  ice  edge  followed  by  seven  parallel  segments  70 
km  from  the  ice  edge.  There  were  great  variations  in  lead 
distributions  between  segments  ranging  from  0.35  to  5.34%  open 
water  (Appendix  B)  with  no  obvious  relationship  to  distance 
from  the  ice  edge.  With  13  data  sets,  73  comparisons  are  pos¬ 
sible,  but  five  segment  pairs  contained  less  than  the  required 
10  elements;  therefore,  only  73  significance  determinations 
were  made.  Sixty-seven  of  these  comparisons  showed  no  signifi¬ 
cant  difference  (P  >  10%).  There  was  a  probable  significance 
level  (P  =  1-5%)  between  segment  pairs  (1  and  12,  3  and  4, 
3  and  12,  4  and  6,  4  and  13)  and  a  statistically  significant 
difference  (P  <  1)  between  segments  1  and  4. 

This  region  is  composed  of  mostly  first-season  ice,  and 
due  to  ocean  swell,  great  fluctuations  in  percentages  of  open 
water  can  occur  rapidly.  Wind  from  the  pack  will  scatter  the 
floes  as  an  opposite  wind  will  compact  the  area. 

EAST  GREENLAND  REGION 

The  East  Greenland  region  includes  the  Greenland  Sea  and 
Denmark  Strait.  The  East  Greenland  drift  stream  represents  the 
major  efflux  zone  of  water,  ice  and  heat  outward  from  the 
central  polar  pack  ice  regions.  This  region  is  often  termed  a 
dynamic  "Ice  Factory"  because  it  is  subject  to  nearly  instan¬ 
taneous  response  to  wide  variations  in  wind  speed  and  direc¬ 
tion.  Great  quantities  of  new,  first-year,  and  sometimes 
multi-year  ice  are  advected  into  the  warm  waters.  In  mid¬ 
winter  months  ( December-March ,  inclusive)  thousands  of  square 
miles  of  new  ice  are  continually  forming. 

The  ref lectometer  data  collected  in  this  region  was 
gathered  during  adverse  weather  conditions  with  respect  to  the 
system  potential.  Therefore,  a  direct  count  of  open  water 
areas  could  not  be  obtained  with  confidence.  Figure  3  shows 
the  percentage  of  ice  concentration  along  the  marginal  ice 
zone  with  respect  to  latitude  for  average  minimum  and  maximum 
conditions.  The  figure  was  constructed  from  data  available 
from  the  Navy-NOAA  Joint  Ice  Center,  Navy  Polar  Oceanography 
Center,  Suitland,  Md .  ,  in  their  "Analysis  of  Eastern  Ice 
Limit,"  1973  through  1980. 
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Figure  3.  Change  in  percentage  of  ice  concentration  with  latitude 
in  the  East  Greenland  Region  during  average  minimum 
and  average  maximum  conditions 


Minimum  conditions  occur  in  the  second/third  week  in 
September.  During  this  time  the  ice  limit  has  retreated  to 
appr oximately  71°N  latitude,  during  which  severe  conditions 
can  leave  a  narrow  band  of  very  low  ice  concentration  along 
the  coast  to  64°.  From  the  southernmost  edge  of  the  ice  limit 
in  the  Greenland  Sea  to  North  Spitzbergen  (approximately 
8 0 °N ) ,  the  data  showed  a  linear  decrease  in  the  percentage  of 
open  water  with  latitude  (Fig.  3). 

At  maximum  ice  extent  in  late  April/early  May,  the  ice 
limit  along  the  east  coast  of  Greenland  extends  beyond  the 
southern  tip  due  to  the  cold  east  Greenland  current.  Data 
analysis  shows  a  high  concentration  of  thin  first  year  ice 
(approximately  70-80%)  from  59°N  to  64°.  The  concentration  in 
this  area  can  change  dramatically  over  a  short  period  of  time, 
particularly  if  the  sea  swell  breaks  the  ice  and  prevailing 
southerly  winds  and  currents  transport  it  away  from  the  coast. 
A  lower  ice  concentration  (60-70%)  is  encountered  from  64°N  to 
72°N  due  in  part  to  the  physical  changes  along  the  Greenland 
coast  and  the  geographic  position  of  Iceland.  From  75°  to  80°N 
a  linear  increase  of  ice  concentration  is  displayed  with  a 
similar  rate  of  change  for  that  area  during  minimum  conditions 
but  with  an  average  15%  higher  concentration. 


WEST  GREENLAND  REGION 


The  West  Greenland  Region  includes  Baffin  Bay,  Davis 
Strait,  and  the  Labrador  Sea.  This  region  parallels  the  East 
Greenland  Region  in  that  it  is  virtually  ice  free  in  the  late 
summer  months  with  high  concentrations  of  thin  ice  after 
freeze-up. 

Analysis  of  ice  edge  movement  in  Baffin  Bay  from  1973 
through  1980  shows  that  minimum  conditions  occur,  on  the  aver¬ 
age,  during  the  first  two  weeks  in  September.  Davis  Strait  and 
Baffin  Bay  become  ice  free  with  occasional  bergs  entering  via 
Kennedy  Channel  or  Lancaster  Sound.  Following  exceptional ly 
cold  winters  or  during  exceptionally  cool  summers  low  ice 
concentrations  will  remain  along  the  coast  of  South  Ellesmere, 
Devon,  and  North  Baffin  Island. 


Maximum  ice  extent  occurs  during  the  last  two  weeks  in 
April.  Due  to  coriolis,  wind,  and  the  cold  Labrador  current, 
ice  growth  proceeds  along  the  western  coastline  of  this  region 
to  a  southern  extent  of  45°N  latitude.  The  ice  near  the  south¬ 
ern  limit  (approximately  45o-50°N)  is  generally  less  than  30 
cm  thick  with  a  rapid  increase  in  concentration  with  latitude 
as  shown  in  Fig.  4 . 
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Figure  4.  Change  in  percentage  of  fee  concentration  with  latitude 
in  the  West  Greenland  Region  during  maximum  conditions 


From  50°  to  the  northern  extent  of  the  ice  limit  at  maximum 
conditions  (approximately  67°)  the  rate  of  increase  of  the  ice 
concentration  from  the  edge  to  100  miles  into  the  pack  decrea¬ 
ses  as  the  average  thickness  increases  ranging  from  30  to  120 
cm. 
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CONCLUSION 


The  use  of  ref lectometer  signal  analysis  for  open  water 
identif ication  in  Arctic  sea  ice  has  proven  successful.  Its 
future  use  regarding  ice  type  identification  with  respect  to 
roughness  shows  great  potential  as  shown  in  Figure  1. 

The  frequency  and  percentage  of  open  water  and  thin  ice 
areas  of  individual  data  sets  and  geographic  regions  are  list¬ 
ed  in  Appendix  B.  The  instability  of  the  'retie  pack,  particu¬ 
larly  in  late  summer  months,  can  lead  to  great  variations  of 
ice  conditions  over  short  time  periods  and  distances  as  shown 
in  the  data.  The  thinner  first-year  ice  of  the  Bering  Sea 
showed  no  obvious  relation  between  distance  from  the  ice  edge 
and  lead  characteristics  with  an  overall  first-year  ice 
concentration  greater  than  98%. 

Within  the  limits  of  this  study  (ice  edge  to  160  km),  the 
Lincoln  Sea  and  North  Slope  regions  also  displayed  no  apparent 
relation  of  open  water  percentages  to  distance  from  ice  edge 
with  ice  concentrations  greater  than  95%.  However,  both 
regions  statistically  yielded  significantly  different  lead 
characteristics  in  data  sets  just  beyond  the  160  km  limit. 
This  may  indicate  a  transition  from  the  marginal  ice  zone  to 
central  pack  ice  with  respect  to  open  water/thin  ice  and 
required  further  investigation. 
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Appendix  A.  Listing  of  Wilcoxson  Sun  of  Ranks  program 
written  by  A.  W.  Lohanick,  NORDA  Code  332 
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APPENDIX  A 


10  *  THIS  PROGRAM  PERFORMS  A  m  WILC  O.'iOM  S  SUN  OF  PfllU  S  TEST  '  AS  DESCRIBED  IN 

"  PRACTICAL  STATISTICS  “  BV  RUSSELL  LANGLEY  pp .  1  68 ,  1 69 ,  DOVER  PUB. 

20  REAL  XI  .  1 :  500  >  f  XI*  1  :500  •  ,Xc  •  1  :  1000  ■  ,  F  1  agl  .  1:500  ■  ,  F  1  492-  1  1*500  -  ,R*nk  <  1  :  1000  ‘ 

30  I NTEGER  P ,  Comb_F  1  1  <_i  1  ,  I r»c  de  •:  •  N  •*,  Nt> ,  Hr  ,  R 

40  DIM  TIM:  1000  til  3  ,  T  4.3 1  $  •  1 :  500  •  C  H  ,  T  4-g2f  *  1  :  500  •  £1  3 ,  Sf  I  30  3 

50  PRINTER  IS  0 

60  1  .  DATA  STATEMENTS  . . . 

70  1  <  MARK  END  OF  EACH  DATA  SET  HITH  9999  ■ 

80  1  EXAMPLE:  DATA  -or,  3, 5, 2, 6,  1,9999 

90  » 

100  »  .  INPUT  DATA  FROM  DATA  STATEMENTS  . 

no  » 

120  FOR  T  i  fn« s  l  TO  2 

130  IF  T  i  rti€ s  1  THEN  READ  R6910r.ll 

140  IF  Time *2  THEN  READ  Pegion2f 

150  FOR  A* 1  TO  500 

160  IF  Tir.it  =2  THEN  Secc*».d_t  if.it 

170  READ  .1  •  A  • 

180  T  49 1 1 • A  •  * ” A " 

190  GOTO  220 

200  S«?c  or.d_t  I  we  :  READ  X2'A- 

210  T  4921* A  >s"B" 

220  IF  \Tkh#»1  •  AND  •  : .  1  ■  A  •  *  999  '•*  •  OR  •Tin.«=2*  AND  • X2 > A > *9999 '  THEN  Off 

230  NEXT  A 

240  0  'f : 1 

250  IF  Tim*=l  THEN  F 1 1 1 ze _ 1 =A- 1 

260  IF  Tum**2  THEN  F  1  1  =  '  re_2=A-  1 

270  NEXT  Time 

280  Con.  t<__i  '  l  e_-  i:<«Fi  ze_l  ♦F  1  I  e_  =  1  :«  i’ 

290  •  . 7.  COMBINE  FILES  . 7 . 

300  FOR  N=1  TO  C  on.b_f  1  )  1  26 

310  IF  N<  =F  1  1  e__=  1  :«  T  THEN  :;c*H'«XDN' 

320  IF  H‘*FiU'’si:€  1  THEN  T  $  ■  N  •  *T  49 1  »  •  N  * 

330  IF  M  F»l«  Ji:ej  THEN  X  c  ■-  N  >  *  X  2 1  N  -  F  1  I  <  _  4  1  z  4  __  1  ' 

340  IF  N  Fil*_$ize_t  THEN  T#< N  .*T4$2f • N-F 1  1 e_i7z«  1- 

350  NEXT  N 

360  «  .  SORT  COMBINED  FILE  . 

370  IF  0$s‘fr  THEN  I  nc  dec  =  1 

380  IF  Of = " D "  THEN  Incde-:=0 

390  CALL  V«c  t  or  sor  ♦  _q •  *  1 ,  T I  ■  *  ■ ,  1  * *o.fc_f  1  1  In.  de  c  • 

400  »  . . .  BEGIN  RANI  ING  .  . 

410  R%nk_uai ue* 1 

420  FOR  P*2  TO  Conb  ♦'  1  le_»i  2f 

4  30  IF  :<c  •  P  'Me  •  P-1  •  THEN  Mo  J  1  « 

440  Tien 

450  T  =  1 

460  FOR  G*P  TO  C  on.t_f  1  1  «  _ »  1  ;e 

4  70  IF  .  0  •=  •;  1  0-  I  ■  Them  T=  T+i  *i.«.  r«2is4ti* 

480  IF  •:  •  0  ■  XcO-l’  THEN  Out 

490  NEXT  0 


500  Out : » 


510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 


Sun,  =  0 

FOR  R= l  TO  T 

Sun.*  *  uf.i  +  R  ani  _•*  at  ue 

R  4  n  1  "  4  I  *.4  e  =  R  4  f‘*  1  1 1  4  1  u  6  ♦  1 

NEXT  P 

FOP  S=P-i  TO  P ♦ T - 2 
R  4nk  1  S  1  x Sun.  T 
NEXT  S 
P  =  F ♦ T -  I 

GOTO  O.e  c  X  w«  or  nd 

No_t  1 e :  1 

R*nl  ■  P-1  ***P4»4  j*alue 
R  4r.1  •  ll  >.J «  *R  ant  ■  4  t  ue  ♦  l 

Chec ^  _♦ or _end:  1 

IF  P  Cof»t  _f  1  1  e_  *  1  ze  THEN  Ne  t 
Part  ■  P  •  *P  4*1 1  _• 1  4 1  ue 
GOTO  Results 
N.?*»_p:  NE  T  P 

1  .  OUTPUT  PENULTS  . 

R.ttul  t  »  .  i 

‘  PRINT  Re  9  1  on  1  f  «  Re  9  1  or.  2 1 

PRINT  ‘A  «  "  :  Re  9ivt.1t*  "l  ~  ”  ;  Re  9 1  *r.2t 
FOP  A* 1  TO  69 

rr  A* 1  THEN  PRINT  C HP# '2 32- 


T 


730  PRINT  USING  “  " 

760  IF  ft«6 9  THEN  PRINT  CHRI*.128> 

770  NEXT  ft 

780  PRINT  "  1  D»t  a  "  fcl  un  1  Tally  1 _ Rank  _ |  ft  r>nl->  |  B  r 

»nk 3  |  " 

790  FOR  ft*l  TO  Coinb_f  i  I  e_s  <  zt 

800  IF  T$ »  ft ' *~B ”  THEN  F 

810  IF  FRrtC  T<  Rank  *  ft  *  »«0  THEN  PRINT  USING  ‘‘ 4X,  PDPI',  12X,  A,  I  IX,  DDPD,  13X,  DDBD" 

;  Kc  »  ft  ' ,  T  f  •  ft  ) ,  R  ink  i  ft  > ,  R «r‘il  ‘  ft  * 

8^0  IF  FPftC  T  •  R  *r.k  •  ft  ■  •  0  Then  PRINT  IJ:ING  •■4. DDt'D*  12:'.,  ft.  I  1  X,  DDBti.  D,  I IX,  DP 

DD .  D  "  ;  .  .'c  *  ft  • ,  T  f  <  ft  • ,  R  ink  i  ft  •  ,  R  an  l  •  ft  • 

8  30  Tot  aI_a*Tot  a!_a*Rarik  •  ft  * 

840  GOTO  He  t  _* 

830  £ :  ‘ 

860  IF  FPflCT*  P»nl  ■  ft  ■  '»0  THEM  FF  ItlT  'J5IMG  ” •»:•:,  P0PD,  1  Is ,  ft ,  1 1 DDDD,  24X,  OPUS" 

;  Xc  •  ft  • ,  T I »  ft  * ,  R  ar.k  <  ft  »  ,  R  ar.k  ■  ft  1 

8T0  IF  FPft>:  T»  Rank  •  ft  0  th£n  FRInT  Ui  I NG  *’  4  f  P  p  P  P ,  1 2  X ,  ft ,  l  1 X ,  D  D  D  D .  D ,  2  2  X ,  D  P 

DO .  D  *’ ; ;  ic  •  ft  • ,  T# '  ft  * ,  R  ar.k  •  ft  - ,  R  ank  *  ft  • 

660  Total  b*Total  b*F  ink  ■  ft  * 

890  N*=.  t  NEXT  ft 

900  “  PRINT  CHP#v l 32 • 

910  FOR  ft=l  TO  80 

920  PRINT  USING  “ 

930  NEXT  ft 

940  PRINT  CHRRC128.' 

950  PRINT  USING  "  7X,  DODD ,  8ft ,  X,  DDPD.  7ft,  9;:*  7ft,  X,  DPDPD.  PP,  3X  .  PDPPD.  DD*‘  J  F  i  1  t_* 

i  z<  I.*'  From  ft ,  " ,  F  i  f  e  <  i  ^  2,"  fro  ft*  F"  ,  “  lot  1 1  i  J  "»  Tot  al_4(  To*  il^b 
960"2^alue:  * 

9 70  Ha«F i )e_si2t_l 

980  Nb»P i U  ii2«  2 


990  R*niNv  T.;.t  al_a.  Tot  a1_b  » 

1000  IF  NOT  ‘‘Na-20-  0P~-  Mb  20  ■  ■  Then  T  at- 1  e _1  ook  up 


1010 

IF 

To.  *l_»- »Tot  *t  t  them  Mr*u» 

1020 

IF 

To'il^b  *Tot  il_i  Them  Nr=NB 

10  30 

Z* • Nr* 

•  ll»«Mt»  1  •  -2*F  •  Jif'p.  Mi- Mt.*'  M  i*Mb*  1 

1  3> 

104O 

IF 

fl$S<r.'5.09  THEM  i  »  *  ■■ .  0 . 2  \  " 

1030 

IF 

>  ABE'  2  •  *3.0®  '  HMIi  'HE:'  Z  '  2.  'S' 

THEN 

Sf  * 

"be  • 

>.>e  e  n 

.2*.  »nd  r. 

1060 

IF 

;  •  =2.5:  hmi.  .hi;  ;•  1 . 

THEN 

SF* 

"be- 

ten 

r.  »nd  5V 

I  (TO 

IF 

'  RES'  Z  •  *1.  *6-  HNti  •  HE':-'  Z  ■  1.64' 

THEN 

s«- 

"  ‘•5* 

.  but 

.  10V 

1080 

IF 

HE? '2'  *1.64  THEM  :»»••  jo\” 

109O 

print  using 

-4H, 2  , dcd,  !■  d t' - ;  ■;  * 

1100 

GO  TO  Pr  t nt 

i 

1110 

f  at- 1  <  lookup:’ 

1120 

CHLL  TiOU'  Mi, Mb, P, if  ' 

1  1  30 

Print  s : ' 

1140 

PRINT 

L  I  Mf  1  •  ,  "  Protit-i  1  M  "  tft  a*  tot 

ft  i i »>p 1 e 

>  c 

imc  f 

rom  i am#  p* 

u  l  at  i  or.  is  "t$t 

1150  PRINT  L INv  ? > 

1 1  so  »:m» 

1170  SUB  Tat.:  <■  INTEGER  n  a.  Nb ,  P ,  SI  1 


1180 

BHTft 

2,8,4, 3,0,0 

1190 

drtr 

2,9,4, 3,0,0 

1200 

DftTft 

2, 10,4, 3,0,0 

1210 

DftTft 

2, 11,4, j, 0,0 

1 220 

DftTft 

2,12,5,4,0,0 

1230 

DftTft 

2,13,5,4,0,0 

1240 

DftTft 

2,  14,6,4,0,0 

1250 

DftTft 

2,15,6,4,0,0 

1260 

DftTft 

2. 16,6,4,0,0 

1270 

DftTft 

2,12,6,5,0.0 

1280 

DftTft 

2,18,2,5,0,0 

1298 

DftTft 

2, 19,2,5, 3,0 

1300 

DftTft 

1310 

DftTft 

3 , 5 , 7 , 6 , 0 , 0 

1328 

DftTft 

3 , 6 . 8 ,  2 , 0 , 0 

1330 

DftTft 

3  .  7 , 8 ,  7 , 0 . 0 

1340 

DftTft 

3,8, 9, 8, 0,0 

1350 

DRTR 

3,9,10,8,6,0 

1360 

DRTR 

3,  10,  10,9,0,0 

1 320 

DRTR 

3,11,11,9,6,0 

1380 

DftTft 

3,  12,11, 10.7,0 

1390 

DRTR 

3,13,12, 10, 7,0 

1400 

DftTft 

3, 14, l 3, 1 t , 7,0 

1410 

DftTft 

3,15,13,11,8,0 

1420 

DRTR 

3,16,14,12,8,0 

1430 

DftTft 

3, 17, 15,12, 8,o- 

1440 

DftTft 

3, 18, 15, 13,8,0 

1450 

DPT  ft 

3,19,16,13,9,0 

« 


1460 

DOT* 

3,20, 17, 14, 9,6 

1470 

DOT* 

4,5. 17.11,0,0 

1480 

DOT  fl 

4,8,13,17,10,0 

1490 

DOTO 

4,7,14.13,10.0 

1500 

DOTO 

4.8,15,14,11,0 

1510 

DOTH 

4,9,18,14,11,0 

1520 

DOT* 

4,10,17,15,12,1 0 

1530 

DOTO 

4,11.18,16,12,10 

1540 

DOTO 

4,12,19,17,13.10 

1550 

DOTO 

4 ,  13,20, IS, 13,  1 1 

1560 

DOTO 

4,14,21,19,14,11 

1570 

DOTO 

4,  15, 22,20,  15.  1  1 

1580 

DOTO 

4, 16,24,21. 15.12 

1590 

DOTO 

4, 17, 25,21 , 16.  12 

1600 

DOTO 

4, IS. 78, 77, 16, 13 

1610 

DOTO 

4, 19,27,23, 17,13 

1620 

DOTO 

4,70,78,74, IS, 13 

1630 

DOTO 

5,5, 19,17,15,0 

1640 

DOTO 

5, 6, 70,  IS,  16, 0 

1650 

DOTO 

5, 7,21 ,20, 1 o , 0 

1660 

DOTO 

5,8,23,21, 17. 15 

1670 

DOTO 

5.9, 24. 22. 18,  16 

1680 

DOTO 

5, 10, 26, 2  3, 19.  1 6 

1690 

DOTO 

5, 11,27,24,20, 17 

1700 

DOTO 

5, 12,28, 26,21, 17 

1710 

DOTO 

5, 13, 30,27,22. 1* 

1720 

DOTO 

5, 14, 31 , 28, 22.  18 

17  30 

DOTO 

5,15, 33,29,23, 19 

1740 

DOTO 

5, 16,34, 30, 24, 20 

1750 

DOTO 

5,17,35, 32,25,20 

1760 

DOTO 

5,  18, 37,  33,26,21 

1770 

DOTO 

5, 19, 38, 34, 27,22 

1  730 

DOTO 

5, 70, 40, 35,78. 77 

1790 

DOTO 

6,6,23, 26, 23,0 

1800 

DOTO 

6. 7,29,27,24,21 

1810 

DOTO 

6,3. 31 , 2 9, 25, 22 

1870 

Dft  T  9 

6,  o,  33,  31 , 26.  2 3 

1830 

DOTH 

6, 10, 35, 37,77,74 

1840 

DOTO 

6,11,37, 34,28,25 

1850 

E0Tm 

6 ,  U,  38,  35,  30, 1 5 

1860 

DOTO 

6  ,  1  3,  40,  37,  3 1 , 26 

1870 

DOTO 

6,  14, 42, 38,  32,27 

1880 

DOTO 

6,  15,44,40, 35,78 

1890 

DOTH 

6 ,  1  6 , 4  6 , 4  2 , 34,29 

1900 

DOTO 

6,  1  7,  4  7,  4  3,  36,  30 

1910 

DOTO 

6, 18,49,45, 37, 31 

1920 

DOT* 

6,19,51, 46, 38, 37 

1930 

DOTO 

6, 70, 53,48,  39.  33 

1940 

DOTO 

7,  7,  39,  36,  37,79 

1950 

DOTO 

7,8,41,33,34, 30 

1960 

DOTO 

7,9,43,40, -5, 31 

1970 

DOTS 

7,10,45,42,37. 33 

1960 

P8T8 

7,  M  ,  4  7,  44,  38.  .  4 

1990 

DOTO 

7, 12,49, 46, 40, 35 

2000 

DOT  9 

7,  1  3, 52, 48, 41  ,  36 

2010 

DOTO 

7 ,  14,54,50,4  3,  37 

2020 

DOTO 

7,  15,  56, 52, 44,  33 

2030 

DOTO 

7,  16, 56,54, 46, 39 

2040 

DOTO 

7,17,61, 56,47,41 

2050 

DOTH 

7,  18, 63,58,49,  42 

2060 

DOTO 

7 ,  1  9 ,  6  5 , 6  0 ,  5  0 ,  4  3 

2070 

DH  T8 

7, 70, 67, 67 , 57. 44 

2080 

DOTO 

8,8,51,49,43, 40 

2090 

DOTO 

8,9,54,51,45,41 

2100 

DOTO 

8, lO, 56, 53, 47. 42 

21  10 

DOTO 

8,11,59,55,49,44 

2120 

DOTO 

3, 12,62,58,51 , 45 

2130 

DOTO 

3, 13,64,60,53,47 

2140 

DOTO 

8,  14,67,62,54, 48 

2150 

DOTO 

8. 15, 69, 65,56, 50 

2160 

DOTO 

8,  16, 72, 67, 58,51 

2170 

DOTS 

8, 17,75, 70, 60, 53 

2  1  80 

DOTO 

•3,  1  77,  72,62,  54 

2190 

DOTO 

3,19,80, 74.64,56 

2200 

DOTO 

2210 

DOTO 

9, 9,66, 62, 56,52 

2220 

I ihTh 

9 ,  I  0 , 6  9 .  6  5 , 5  S ,  5  8 

2230 

DOTO 

9, 1 1 . 72,68,61 . 55 
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DATA  9,  12, 75, 7  1 ,oJ,b. 

DATA  3,  1  3,  78,  7  3, 65,  59 
DATA  9, 14,81 , 78,67,88 
DATA  9, 15,34,79,69,82 
HATH  9, 16, 87, 82, 72, £4 
HATH  9,  17 ,90, 64, 74,6c- 
DATA  9, 18,93,87 , 76,63 
DhTA  9,  19.  96, 90,  78.  7£i 

DATA  9, 20. 99, 93. 81 , '1 
DR  V  9  10,10,82.78,71. 65 
DATA  10, 1 1 , 86, 81 , 73, 67 
DRTR  10, 12,69,84,76.69 
DRTR  10,13,32,88,79.72 
DRTR  10, 14,96,91,81, 74 
DRTR  10,15,99.94,84,76 
DRTR  1 0, 16, 103. ?7, 36 , 78 
DRTh  10,17,106,100.89,80 
DRTR  10,18,110,103.92.82 
DRTR  10,19,113,107,94, 84 
DRTR  10,20,117,110.47,87 
DRTR  11,11,100.96,87,81 
DRTR  11,12, 104,99.90,83 
DRTR  1  1 ,  1  3, 1G8,  103. 93,86 
DRTR  11,14,112,106,96,68 
DRTR  11,15,116,110,99.90 
DRTR  1  1  ,  16,  120.  1  1  3,  102, 93 
DRTR  1  1.17.  123.  1  17.  105,  95 
DRTh  1 1, 18, 127 , 121, 108,9* 
DRTR  11,19,131,124,111.100 
DRTR  1 1 , 20, 1 35, 128, 1 1 4, 103 
DRTh  12,12,120,115. 105.9? 
DRTR  12, 13, 125, 119, 109, 101 

DRTR  12. 14, 129, 123, 1 12, I03 


15, 

:o, 

127, 

115, 

10* 

1  *, 

133, 

131. 

117, 

109 

17, 

142, 

135, 

122, 

1  12 

18, 

14*, 

139, 

125, 

115 

19, 

150, 

143, 

129, 

118 

20, 

155, 

147, 

132, 

120 

1  3 , 

142, 

136, 

125, 

1  1  7 

14, 

147, 

141  , 

129, 

120 

15, 

152. 

145, 

13  4, 

12  3 

16, 

15*, 

150, 

1  3* , 

12* 

17, 

1*1  , 

154, 

140, 

129 

13, 

1  ** , 

158. 

144, 

1  3  3 

19, 

171  , 

1*3, 

148, 

13* 

20, 

175, 

1*7, 

151  , 

1  39 

l**, 

1*0, 

1  47, 

1  37 

1% 

in. 

1*4  , 

l  5  1  , 

1  4  I 

16, 

17*, 

1*9, 

155, 

t  44 

17, 

182, 

174, 

159, 

1  48 

1  8 , 

137, 

179, 

1*3, 

151 

19, 

192, 

183 

1  *3 

155 

20, 

197 

133 

172 

159 

15, 

192 

184 

171 

1*0 

16 

1  97 

190 

1  75 

1  *3 

,  17 

20  3 

195 

1  80 

1  *7 

,  18 

62  0  8 

200 

134 

1  7  1 

.  19 

214 

205 

1  89 

175 

,  20 

220 

,210 

1  9  3 

1  79 

.  1  * 

219 

21  1 

19* 

134 

,  1  7 

225 

217 

201 

1  38 

,  is 

231 

,  i-  2  — ’ 

20* 

1*2 

.  19 

2  37 

228 

210 

19* 

,  20 

24  3 

.234 

215 

201 

.  1  7 

249 

,  240 

22  3 

.210 

,  19 

255 

,  24* 

2  2  8 

,214 

.  19 

2*2 

,252 

2  34 

,219 

,20 

2*3 

,  258 

2  3  9 

.  22  3 

.  18 

280 

,270 

252 

,  2  -■  < 

.  19 

23  7 

,  277 

6.58 

,  242 

,20 

294 

,  2  8  - 

2*  3 

.  24  7 

.  19 

213 

,  303 

23  3 

,  2*7 

•  20 

326 

,  309 

289 

,2  72 

DRTR  20,20, 348. 337, 315, 298 
FOP  fi-1  TO  182 

READ  HI ,N2,P1 .P2.P3.P4 

IF  NOT  (iNl'Ill  ■  AND  ■  Hb=H2 


OP  ' Ma*M2  >  AND  Nb-Hl 


t_4 

3020 

30  3O 

3040 

3050 

306C 

3070 

3060 

30*0 

3100 

31  10 
3120 
3130 
3140 
3150 
3  1  60 

3iro 

3160 

31*0 

3200 

3210 

3220 

32  30 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 

35  30 
3540 
3550 
3560 
3570 
3530 
3590 
3600 
36 1  0 
3620 

36  30 
3640 
3650 
3660 
3670 
36S0 
3690 
3700 
3710 
3720 
3738 
3740 
3750 
3760 
3770 
3780 


IF  R  F4  THEN  il  =  "  0.2V 

IF  >  R  =  P  4  MNP  *  R- F  •  ThEN  Sl  =  "t*«  tween  .2’.  and  1*. 

IF  R  =F  1  ■  HMD  ■  R  F21  ThEN  51=  “b«  t  wee  n  \%  and  5’." 

IF  •  R  *P2  ■  HMD  F  F  1  •  THEN  £.l  =  ,‘  5*.  bu»  <  10V 

IF  R  =P1  THEN  if*"  10V 
SUFE"  I  T 

Ne  HE  T  H 

6*  *  "  Not  in  t  it  I  e  " 

iUPENp 

iUE  v«  c  t  or  £  or  *  _g  ‘  H  ■  *  * ,  Af  •  *  • .  INTEGER  I  1  ,  J  1  ,  I  nc  de  C  • 

INTEGER  Logt  uo 

n* n ♦ i - 1  i 

Logt  wo*  I  NT  «  UGT  v  N  .•  LGT-2  +1 

CALL  Qsort  <  A<  ♦  Af  <  ♦  >,Log»  ■•«*,  11.  T  3  ,  Inc  dec  • 

>UE>E>  I T 


>UF  0  5  or  t  INTEGER  Log,  II,  J  t  •  Inc  dec  • 

OPTION  BASE  1 

JIN  L ^ Log  * , U* Log  * 


N*  1  • 

1  =  11  1 

J=J1  • 

jt  art  1;  IF  J  a  j  THEN  Ne  f  gr  •.■up 
St  ar  »  2:k  =  I 

1 2=  I  NT  •  1  «  2  » 

T  =  Fk  12  ' 

Tf=At<  12  .•  •  ' 

IF  Inc  dec =0  THEN  PI 
11 J  IF  A  <  l  >  =  T  THEN  Lowmiddlel 
GOTO  3320 

31 :  IF  A  <  I  > '  =  T  THEN  L o »♦. i  d d  I  e  1 

AC  12 >=*<!>  ' 

A ».  J  >  *  T  ' 

T=ftU2)  1 

A  $  •  1 2  >  =  A  *  • -I  •  • 

AfvI..»T*  ! 

T  f  =  A I  v I2>  ' 

Lour,,  i  dd  1  e  1 :  L=J 

IF  Incd€-:=0  Then  P2 

12:  IF  h>  J  .  =T  THEn  Middlehigh 

GOTO  3430 

>2:  IF  A >:  J  1  =  T  THEN  Middlehigh 

A<  1 2  >  *  A  ■  J  • 

A  <  J  >  =  T 
T  = A  1  12  • 

Af < I2’=Af •  J  • 


:  e  *  stack  po  nit  ir  , 
Get  l  o  er  e  n  d  i:«  i  n  t . 
Se  ♦  uppe  r  e  ndpo  i  n  t  . 


Determine  the  midpoint  of  a  segment. 


(heel  to  i e e  it'  1  owe r  e ndp o i n t  an d 
midpoint  are  in  order.  If  not, 
switch  them. 

Reset  mi dpo i nt . 


•«’  upper  endpoint. 


1  Che  cl  to  iee  it  he  midpoint  and 
1  the  upper  endpoint  are  in  order. 
1  If  not.  -witch  them. 


A  $  v  J  >  *  T  t  >  — 

T$=AI t  12  ^  ' 

IF  I n c  d e ■:  =  0  THEN  p  ?• 


t  3: 

IF  A r I  .  =  T 
GOTO  35  j0 

THEN 

N i dd 1  eh i gh 

3  3 : 

IF  A I  >  >  =  T 

THEN 

M  i  dd  Uhi  gh 

A 1  12  >**•  I  ■ 

A'  I  ■  =T 
T-ft t 12  • 

At •  I2'^=Ht'  I  • 

fit •  I  '  =T t  » 

Tt  =  At»  12  ■  1 

Middle  high:  L  =L -  1 

IF  I n c d e c  =  0  Then  p 4 
f  4  :  I F  ft <  L  /  ■  T  THEN  M  i  dd  f  e h  »  gh 

GOTO  3640 

34:  IF  A-L-  T  THEN  Middlehigh 

T  1  =  A  ■  L  ■ 

T 1 1  =  Ht ■ L  *  i 


'  Check  to  see  if  'he  s w i t c h  1 e f t 
'  the  lower  endpoint  and  the  mid- 
1  point  m  order. 

*  If  not,  switch  them. 


1  Decrement  the  upper  endpoint. 


1  Check  t iee  if  *he  new  upper 
*  endpoint  j£  >n  order. 


?t  ►.pup: 


15: 

35: 


\  =\.  ♦!  I 

IF  I ncdec  =0  THEN  P5 
IF  A- h > - T  THEN  it e pup  ' 
GOTO  3710 

IF  A • K ) > T  THEN  St e pup 
IF  K ' L  THEN  Passed  * 

A  v  L  >  *  A  ■  k  •  ' 

A  <  K  )  *  T  J  » 

Af  <L  5  =A*  '.K  >  » .*.  • 


If  not,  i  a-  e  the  upper  endpoint  and 
increment  the  lower  endpoint.  Now 

check  if  the  lower  endpoint  is  less 
t  h an  t  he  m i dpo inf .  IF  not ,  t hen  switch 
t  he  upp e  r  and  1 o w e  r  endpoints. 


At  t  K  P  =  T  1 J  ‘  /  r 

GOTO  Middlehigh 

Passed:  IF  L  - 1  =  J  -  k  THEN  Sforehigh  1  Sort  the  shortest  segment  firs*  . 
L i M  >  *  I  1  St  or e  t  he  1  owe r 


18 


3790 

U<f1)  *1 

1  endpoint*. 

3800 

I-K 

1  Set  the  new  lower  endpoint. 

3810 

n»nM 

1  Push  the  stack 

3820 

GOTO  3870 

38  30 

itorehigh!  LaH'*K 

1  Store  the  upper 

3840 

LK  M  >  *  J 

1  endpoints 

38*0 

J*L 

1  Set  the  new  upper  endpoint. 

3880 

rt  *  m  ♦  i 

•  Push  the  stack 

3870 

IF  J-I  =11  T  ME  M 

St  art  2 

3880 

IF  I  =  1 1  THEN  it 

art  1 

3890 

3900 

Inc:  I-I-l 

1  Increment  lower  endpoint. 

3910 

IF  I»J  THEN  tU  t group 

'  It  the  current  segment  >*  sorted,  then 

3920 

T*fl* I*t  > 

3830 

T$*ftf *  l*\  •  1 

3940 

IF  I  nc  d«  C  =0  THEN  lie- 

3350 

16*.  IF  ft'  I  >  *T  THEN  Inc 

3960 

GOTO  3980 

3970 

j6:  IF  ft< I 3  >*T  THEN  Inc 

1  Check  to  see  if  ne-xt  element  is  in  order 

3960 

K«  I 

1  Insert  element  in  otherwise  sorted  list. 

3990 

Copy!  ftO. ♦l>*ft,>  * 

'  This  section  dumps  the  array*  up. 

4000 

fti  •  y  1  =  ft  f  •.  k  »  ' 

.  ✓  ' 

4010 

K*k- 1 

1  Prepare  to  bump  ne  t  element. 

4020 

IF  I nc d€ C “0  THEN  D 7 

4030 

17:  IF  Tc  ft  <  k  >  THEN  Cop' • 

4040 

GOTO  4060 

4050 

37 :  IF  T  ftN.  •  THEN  Cop-- 

1  C  h  e  •;  t  »  o  if  insert  ion  is  here. 

4080 

ft <  h ♦  l  ' ®T 

•  It  so.  f  he n  insert. 

4O70 

ft  t  •  k  *  1  *•  *  T 1  < 

408O 

GOTO  Inc 

4090 

He  t group:  M*M- 1 

1  Pop  tt.<  ;i  act  . 

4  1O0 

IF  M  =  0  THEN  Out 

1  C  he  c 1  for  e  nd  c  end 1 1  ion*. 

4110 

I  *L  1  M  ' 

1  Pe  s  t  ot  e  t  he 

4120 

J*U*  M 

1  pre*  lOuf  endpoints. 

4130 

GOTO  3870 

4140 

Out:  SUBEXIT 
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APPENDIX  B.  Class  interval  in  meters 


APPENDIX  B.  Class  interval  in  meters 


23 


15 .862  I  2375 


APPENDIX  B.  Class  interval  in  meters 
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Appendix  C 


Regional  and  overall  histograms  showing  lead 


width  and  frequency  distributions 


H'  1111:1  >.'(  Um  I  i 


APPENDIX  C 


Lead  width  and  frequency  distribution  from  figure  2  track  5 


O  O  O  O  O  O  C-  G  C  G  G1  - 
'■»*.  ll'  ,rr'  -''j  u"  'i'.'  ^  t f  •  ■  j  c  c* 


m  O  Q  I  i  f  f  i  »  i  < 

>*’*>  'X'  Q  O  O  Cj  O  O  O  G 
Cr'  «*'j  u“  C  j  -  t  r  .  G> 

~  ~  ~  Oj  j  .-.J  .V, 

Lead  width  and  frequency  distribution  from  figure  2  track  1 
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APPENDIX  C 


4  0 


i-  T  h  .  .t; 


Lead  width  and  frequency  distribution  from  figure  2  track  2,3,4 


t 

,:LL!:n 

’•--‘-i-  i 

1  e  D  E:  j 

1 4 oL'  ! 

t. 

!  2  0  L  , 
ico[:  i 

6  0  i 

i 

_  _  ►  I  1— 

cfUL'  i 

JJ^ 


0,'EFRLL 
C.!  =  T  =  : BUT! on 


O  O  O  O  O  O  Ci  O  O  O  O 
'jD  <?"•'•!  u  10  — •  *T  {  -  Z  —  O 

[  i  (  —  — ■  ■ —  •  •  J  ■  •  J  •'•J  '<£• 

in  oo  1  1  ’  •  1  1 
— •  ro  ul1 0  O  O  O  O  O  C1 0 
ij>  .--j  u*.  -  -r  r .  c 

—  — ■  —  r-j  .--j  .*-j 


Lead  width  and  frequency  distribution  from  figure  2,  all  track  lines  combined 


27 


Appendix  D.  Example  of  HP-9845B  printout  from  Wilcoxson 


Sum  of  Ranks  program  listed  in  Appendix  A 


PRECEDING  £*0®  BLANK -NOT  TlLim) 
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APPENDIX  D 


fl  *  Nor  Slo  Oa  B  »  Nor  Slo  Ob 


1  Tally  I 

Rat'il;  *  *  a  1  u€  £ 

1  H  f  at'i  l  =  | 

B  r  aril  = 

63 

ft 

i 

i 

72 

B 

2 

2 

86 

B 

3 

3 

96 

B 

4.5 

4.5 

96 

ft 

4.5 

4.5 

120 

B 

€ 

216 

ft 

>■ 

*7 

( 

480 

ft 

8 

8 

528 

fl 

9 

9 

624 

ft 

10 

10 

■’68 

fl 

1  1 

1  1 

864 

ft 

12 

12 

1392 

ft 

13 

13 

1440 

ft 

14 

1  4 

1’28 

B 

15 

15 

4080 

fl 

1  6 

16 

1 1  f  r  o  m 

8, 

5  from  B 

Tot  al  £■: 

105.50 

30. 50 

Pr ob ab i 1 i t  y 

t,  h  at 

hot.  h  s  am y.< '  >5  s 

c  am  €  t  r  o  m  s  a  m  6 

popti  1  at  i  on 

is  >10 

Nor  S 1 0  6a 


B  =  Nor  Slo  Ob 


Bata  values- 

1  Tally  I 

[  h  rank*- 

|  B  rank: 

408O 

ft 

l 

i 

r’23 

E 

•y 

2 

1440 

ft 

o 

3 

1392 

ft 

4 

4 

364 

ft 

5 

5 

76S 

fl 

6 

6 

624 

ft 

7 

7 

328 

ft 

o 

8 

480 

fl 

9 

9 

216 

ft 

10 

10 

120 

B 

1 1 

1 1 

96 

ft 

12.  5 

12.5 

96 

B 

12.5 

12.5 

86 

B 

14 

14 

72 

B 

15 

15 

63 

ft 

16 

16 

11  from  fl,  5  from  B  Totals:  81.50  54. 50 

Probability  that  both  samples  came  from  same  population  is  >  1 0 
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